MA MA 


第 28 卷 第 3 期 工 程 数 学 学 报 Vol. 28 No. 3 
2011 年 06 月 CHINESE JOURNAL OF ENGINEERING MATHEMATICS June 2011 








文章 编号 :1005-3085(2011)03-0307-08 
求解 鞍点 问题 的 多 项 式 加 速 超 松 弛 方法 * 


WEF, FAEN? 
(1- 浙江 工业 职业 技术 学 院 ， 绍 兴 312000; 2- 华东 师范 大 学 数学 系 ， 上 海 200062) 
摘 要 : 为 了 快速 有 效 地 求解 大 型 稀疏 鞍点 问题 ， 在 广义 逐次 超 松弛 (GSOR) 迭代 算法 的 基础 上 ， 结 
合 Chebyshev 多 项 式 加 速 技术 ， 本 文 构 造 了 一 种 多 项 式 加 速 超 松弛 迭代 算法 ， 并 研究 了 该 算法 
的 收敛 性 ， 通过 讨论 加 速 后 迭代 矩阵 的 收敛 性 证 明了 新 方法 比 加速 前 的 迭代 法 具有 快 的 收敛 速 
EE. 数值 例子 也 表明 新 方法 提高 了 GSOR 算法 的 收敛 效率 . 
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分 类 号 : AMS(2000) 65F10; 65F50 中 图 分 类 号 : 0241.6 文献 标识 码 : A 
1 引言 


大 型 稀疏 鞍点 问题 的 求解 在 很 多 应 用 领域 都 出 现 过 ， 并 且 非 常 重要 ， 如 计算 流体 力学 中 
的 Stokes 方程 和 电磁 学 Maxwell 方 程 的 有 限 元 离散 以 及 二 阶 椭圆 方程 问题 的 混合 有 限 元 方法 ， 
带 有 限制 条 件 的 二 次 优化 问题 ， 线 性 弹性 力学 问题 ， 图 像 处 理 问题 ， 最 小 二 乘 问题 等 中。 本 
文 主要 考虑 具有 如 下 形式 的 鞍点 问题 


A a 


其 中 4 € Rmxm 为 对 称 正定 和 矩阵，B € R” 为 列 满 秩 和 矩阵 ， 即 rank(B) =n, HEr, pE 
Rm， 向 量 y, 9 e R^ (这 里 我 们 假定 m > n), BT Æ BHH, p, g 是 已 知 向 量 . 在 这 
样 的 假设 下 ， 鞍 点 问题 (1) 有 唯一 解 由 ， 当 前 已 经 存在 很 多 方法 求解 鞍点 问题 ， 其 中 包括 
直接 法 、Krylov 子 空间 方法 、Uzawa 算 法 、SOR-like 算 法 四 以 及 它 的 推广 算法 个、HSS 算 
法 、GAOR 算法 趾 等 .由 于 直接 法 在 求解 大 型 稀疏 线性 方程 组 时 ， 可 能 会 产生 许多 非 零 元 ， 从 
而 导致 运算 量 大 量 增加 ， 而 且 对 于 条 件数 很 大 的 系数 矩阵， 直接 法 往往 会 产生 很 大 的 误差 . BM 
此 在 求解 大 型 稀 朴 线性 方程 组 时 ， 我 们 通常 会 选择 迭代 法 . 迭代 解法 简单 有 效 ， 只 需 较 小 的 存 
储 空 间 ， 容 易 实现 ， 从 而 在 解 鞍点 问题 的 一 类 方法 中 迭代 解法 得 到 了 大 量 的 研究 结果 . 

鞍点 问题 (1) 是 对 称 不 定 的 、 系 数 矩 阵 对 角 块 中 含有 奇异 矩阵 的 线性 系统 它 的 这 些 特 点 使 
得 一 些 著名 的 迭代 算法 ， 如 PCG AEM SOR HAS 难以 直接 应 用 ， 因 此 ， 有 许多 学 者 提出 了 
一 些 新 的 迭代 算法 用 于 求解 这 类 问题 . 其 中 Uzawa $1109! 和 预 条 件 Uzawa 算法 Lt9 是 一 类 较 早 
被 提出 用 于 求解 鞍点 问题 的 经 典 欠 代 算法 . 经 典 的 SOR 迭代 是 非常 简单 而 有 效 的 方法 之 一 ， 
在 工程 计算 中 有 着 广泛 的 应 用 . 鞍点 问题 (1) 的 特点 使 SOR. 算 法 B9 难以 直接 应 用 ，Golub 等 
人 通过 对 系数 矩阵 进行 适当 的 分 解 ， 构 造 出 了 一 种 SOR-like 和 迭代 方法 四， 在 此 基础 上 白 中 治 等 
人 又 提出 了 一 种 修正 后 的 GSOR 迭代 算法 B3 


收 稿 日 期 : 2009-10-09. 作者 简介 : 潘 春 平 (1978 年 3 月 生 )， 男 ， 讲 师 . 研究 方向 ， 数 值 代数 . 
* 基 金 项 目 : 浙江 工业 职业 技术 学 院 科技 计划 (KY2010109). 


308 I E 数 学 学 R 第 28 卷 








为 了 表示 方便 ， 我 们 将 方程 组 (1) 等 价 表示 为 


Cu )6)-U) : 


将 方程 组 (2) 的 系数 矩阵 .4 作 如 下 分 裂 


A B 
Az —-D-L-U, 
-BT 0 


A 0 0 0 0 -B 
D- s ug S,coU ; 


Qc R” 是 对 称 非 奇异 矩阵 设 


其 中 


wM r AANER, Im E€ g"*", In e R*X" 分 别 是 m 阶 和 nn 阶 单位 阵 ， 考 虚 如 下 一 般 
的 SOR 迭代 格式 


pere -1 a) -1 p 3 
em - (D - QL) [1 - Q)D 4 QU] P" 40D:0L)90 NI (3) 


简 记 为 


其 中 


e lA 0 
Hn 5 07 (D - 9L) = 


B? 1 


T 


(1— wd —wA-!B 

Mw) = (D- QL)"! [1 - Q)D 4 QU] = . (4) 
(1—w)rQ-1BT TI 一 wrQ-LB74-1B 

由 此 得 到 GSOR 算 法 : Qe 名 xn 是 对 称 非 奇异 矩阵 ， 给 定 初始 向 量 z(0) c R”, y9 < 


Rn 和 两 个 松弛 参数 w, 7 > 0， 当天 = 0,1,2,….， 直 到 和 迭代 序列 (207,09 yry esi, ir 


| g(k+1) Z (1 ww)z(p) +wA-i(p = By(9)), (5) 


yk+1) 所 一 y? 十 TQ-L(BTz(k+D 一 q). 
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2 M Chebyshev 多 项 式 加 速 的 GSOR 迭代 算法 
假设 入 是 迭代 矩阵 Mon 的 一 个 非 零 特 征 值 ，(z7, j7)7 e R^ 是 其 对 应 的 特征 向 量 ， 则 


由 (4) 可 得 
(1—wA) -wB 元 4 0 
PO s 
0 Q Ü -rBT Q 


1 -w - A)A$ = wBjj, 
? w - JA = uBij (5 


iH 


e 


即 


(A - 1)Qj = MBTE. 
由 于 和 4 是 对 称 正定 的 ， 由 (7) 可 得 (1 -w -Až = wh-1BY， 若 入 关 1 一 w， 由 于 Q@ 是 非 奇 异 
的 ， 进 一 步 可 得 到 (1 一 一 和 (和 一 1)= rwg, 这 里 =Q-1B7h4-1B. A—-1—w£0, 
WE BJ = 0， 一 wQy = ArBTz, uffág = 0, z e N(BT), N(BT) 是 矩阵 B7 的 零 空 间 ， 因 
此 ， 入 = 1 一 w 是 M7) 对 应 于 特征 向 量 (j7,0)7 的 一 个 特征 值 . 当 入 关 1 一 ww 时 , 设 4 是 J 的 
一 个 特征 值 ， 可 得 (入 一 1)(1 一 w 一 入 ) = 和 Arwx， 即 和 满足 二 次 方程 


A? + (rw 二 w 一 2) 入 +1 一 ww=0. (8) 


通过 解 (8)， 我 们 可 得 当 w E [ihr +o) i Meo 的 特征 信和 是 实数 ， 由 GSOR GERA 
法 的 收敛 性 知 ， 当 0 <w <2 且 0 <+ < A9) 时 ， GSOR .迭代 算法 收敛 13 则 当 0 < w < 
2H0«r«20-m, fpei, TAE 


Tey <w<2, 0<T< EX 
» M," 的 特征 值 和 是 实数 且 | 和 | < 1. 这 样 ， 我 们 给 出 用 Chebyshev 多 项 式 加 速 迁 代 法 的 前 
是 ， 记 
e ( a ) x p ) 
则 GSOR 迭代 格式 简 记 为 
utd) — Mou + c. (9) 
18. EXE AU SUM E TARIME, ARRU TF 
ue) = oe [o (Mo, pu +e) + (1 v)u(9] 4- (1 — g*9yu7D, (10) 


其 中 m0) 是 任意 初始 向 量 ， 且 


其 中 
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o, B 分 别 是 和 矩阵 Mi) 的 最 小 和 最 大 特征 值 的 估计 值 ， 由 此 得 到 如 下 Chebyshev 多 项 式 加 
速 的 GSOR 和 迭代 算法 ( 简 记 为 C-GSOR 和 迭代 算法 ): RQ e R"x* 是 对 称 非 奇 异 和 矩阵 ， 给 
初始 向 量 z(0) e R”, y% e R"* 和 两 个 松弛 参数 w, r > 0， 令 k = 12,3,…， 直 到 迭代 序 
F (29^, 09^ )T) 收敛， 计算 

z 0) = (1— w)z (9 -- 5A-1(p - By®), 

y(O = (9 +7Q 1 (BT QO — q), 

g (1) = gU *D[(1 — vu)a( 9 + vo A71 (p - By(9)] + (1 — g**9)g(-9, 

y(F*) = ptD [VD + vrQ-1 (BT (D — g)] + (1 — gy (70, 


3 C-GSOR 和 迭代 算法 的 收敛 性 分 析 


定理 1 设 4eRmxm,Qe R” 是 对 称 正定 矩阵 , B € R"*" 为 列 满 秩 矩阵 , 令 Lin: Hmax 
分 别 为 了 = @-1B7A-1B 的 最 小 和 最 大 特征 值 。 若 参数 ww, 7 满足 以 下 条 件 : 








1) 0<7< LLL M, w-(T) <w <2; 
1 2(2 一 
2) r 55 T 时 ，w+(r)<w< 2; 
m 4T Umi ATL. 
x (r) S (1 RE TH min)" : w«(r) E (1 T fi^ 
证 明 因为 当 dud) 
TL —w 
xs 4^ SUI e 
FI, Mior) 的 特征 值 入 是 实数 且 | 和 | < 1， 又 因为 
A ATHmin —— 
(1) Mücrs LL LM, 对 任意 的 几 ， XS ; 
MAT np 1 2(2—w) w) AT max 
(à) - y U min Umax DoR WhHmax Ris 对 任意 的 did 有 Tr d T E 2 (1 十 Tax)2 
所 以 上 述 定理 成 立 . 


定理 2 w, 7 满足 定理 1 的 条 件 ， 则 C-GSOR 迭代 算法 收敛 . 
证 明 设 


是 算法 的 精确 解 . 令 5 = u —u0, Enc = qc(M(w,n)C 中 .于 是 要 证 明 所 给 出 的 
迭代 格式 收敛 ， 只 需 证 明 p(gs(M(w,7))) < 1 即 可 . H Chebyshev £ W iE MEA g(r) ER 
于 z 的 有 次 多 项 式 . 定义 如 下 


wo «n (5 5)/n (5? 
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其 中 ZL(z) Ek RERA, EXA 
T(z) = 人 ((z+ Vr? equ (z — 2? 1), 一 co < z < +00, 


X Bl X o, 满足 定理 1 的 条 件 时 ，Mw .的 特征 值 是 实数 且 || < 1， 令 和 (M1) € 
læ, 8] C (—1, 1). 注意 到 





27z 一 0 一 acN\| 
ax, [n( 75777 )] 71 (12) 
4£-2 az 2>], 当 z > 1 时 ， Ty (x) )=#((z+ Vz? — 1 ) 十 一 VZ22 — 1)*), 可 得 
£ 2g* 
P(g (Mos7)) max — lau) < mex la) = (T6) = prx <l 09 


其 中 o 二 &2 一 VEIT= (E+ VE -1)9 <1， 所 以 该 迭代 法 在 w,7 满 足以 上 条 件 时 收敛 . 
下 面 给 出 GSOR 迭代 算法 最 优 参数 的 相关 结论 D3， 
” D 
VES n 


vVl1-w, Zi 0<w<w-(T), 
M(2 — o — ropas. toO — Tata — AE o), 
p(M,,r) 一 若 w-(r) sws wo(r), (14) 


$ (Tumax +w — 2 + y(TwHmax +w — 2)? — 4(1— w) ), 


车 wo(r) <w< 2; 


1 2(2 — w) 
jf wm e ue A 


vV1-w, 4 0<w<w4(T), 
p(Mwr) —- 4 l(renumax +w — 2+ / (Tulmax +w — 2)? — 4(1 — 0), (15) 


Xi w+(r)<w<2， 


这 里 
ÅT limin 4T Umax 4 
w-(T) = s, wlt) = ~, wo(T) = ~, 
Sea Oa 100) 7 aa d aac) 2 
因此 ， 最 优 参数 为 
4VHminHmax 1 V Hmax 7 y Hmin 
Wopt = M, ) = 


TU m— AD) T, x REL opt»Topt 
(V Hmin 十 y Hmax y s y Hmin max d iis: V Hmin + Vmax 


定理 3 假设 wopt 和 7opt 为 GSOR 和 迭代 算法 的 最 优 参数 ，C-GSOR. 迭 代 算 法 的 迭代 矩阵 
为 ae(Mw nm)，p(Me no) 为 最 优 谱 半 径 ， 则 


p(ax (Ms isa) ud p (Miu. ,Topt) ) : 
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证 明 Wo-p(Muuu» a= -o =p, WE = 号 —1. 由 


5 








H ^l 
可 得 = iest) 02 
i o Pg Hi 
又 由 于 
p(qk( Mn)) = NE RENE lax (wu) 
< max Jala) = (n0) ^ = og <1 (16) 
一 a<z<8 k 14 c?k 3 
. audios di. 
p (Musso) lu p TL (- + r) mE (1 E 02)R， 0 <o< T5 (17) 
用 数学 归纳 法 可 证 
20 2^g 
EL GM" ERE Ao MEOS 
LoF < (gne 0coc«1, (18) 
即 
p(ak Mu. sus) < pr Magus) (19) 
4 数值 例子 


下 面 给 出 应 用 本 文 算法 的 数值 例子 . 设 4, BB 分 别 是 一 个 m x m, mm x n 的 矩阵 ， 满 足 
i+1, i=j, 

A= (Qi,j )mxm = 1, i-—3j|l21, B= (bi ;)mxn = | 
0 ”其 它 . 


HE p, q 可 以 任意 选择 . 对 于 SOR-like 方 法 ，GSOR 方 法，C-GSOR 方 法 ， 我 们 均 选择 预 处 
HAERE Q = BTB. 计算 终止 的 条 件 是 


|lrxll> _ 106 ux R e g 
l[rolla Sac —g | | -BT 0 ) | y) d (20) 


这 里 SOR-like 方 法 和 GSOR 方 法 的 最 优 参 数 根据 已 有 的 结论 确定 . 下面 给 出 数值 实验 的 结 
果 ， 所 有 的 结果 都 是 在 Intel E2180 2.0GHZ CPU, 2.0G Memory 条 件 下 获得 的 . 

从 表 1 至 表 5 中 可 以 看 出 ， 在 同一 精度 的 条 件 下 ， 将 加 速 后 的 新 算法 C-GSOR 与 加 速 前 
的 GSOR 算法 和 SOR-like 算 法 相 比 ， 当 w, 7 同时 取 最 优 参 数值 时 ， 新 算法 C-GSOR 对 GSOR 
算法 有 一 点 加 速 ， 但 加 速 的 效果 不 是 非常 明显 ， 新 算法 C-GSOR 比 SOR-like 算 法 具有 更 快 
的 收敛 速度 和 更 少 的 迭代 次 数 ， 当 w, 7 不 选择 最 优 参 数值 时 ， 新 算法 C-GSOR 能 较 有 效 的 
加 速 GSOR 算法 和 SOR-like 算 法 收敛 速度 . 总 的 来 说 ， 新 算法 C-GSOR 能 够 对 GSOR 算法 
和 SOR-like 算 法 有 一 定 的 加 速 ， 新 算法 提高 了 GSOR 算法 和 SOR-like 算 法 的 收敛 效率 ， 特 别 
是 对 于 SOR-like 算法 具有 一 定 的 竞争 优势 . 


j i-jtm-n, 
0， 其 它 . 
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表 1: SOR-like 算法 收敛 迁 代 次 数 和 时 间 


mu [ ewe, 
150 1.9533 1032 











A2: GSOR .算法 收敛 迭代 次 数 和 时 间 


CPU (4) 





表 4: GSOR 算法 收敛 和 迭代 次 数 和 时 间 






表 5: C-GSOR 算法 收敛 迭代 次 数 和 时 间 


CPU ( 秒 ) 
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The Polynomial Acceleration of GSOR Iterative Algorithms for Large 
Sparse Saddle Point Problems 


PAN Chun-ping!, WANG Hong-yu!? 


(1- Zhejiang Industry Polytechnic College, Shaoxing 312000; 
2- Department of Mathematics, East China Normal University, Shanghai 200062) 


Abstract: In order to solve large sparse saddle point problems (SPP) quickly and efficiently, we con- 
struct in this paper a polynomial accelerative iterative algorithm through accelerating the generalized 
SOR (GSOR) iterative algorithm and using the Chebyshev polynomial. The convergence of the algo- 
rithm is also studied. Meanwhile, by examining the convergence of accelerated iterative matrix, the 
result shows that the new method converges faster than the GSOR method. Numerical results further 
demonstrate that the new method is more efficient than the GSOR method. 

Keywords: saddle-point problems; iterative method; GSOR method; Chebyshev polynomial 
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